Conductance quantization and snake states in graphene magnetic waveguides 
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FIG. 1: 


Magnetic field profile ((2} for magnetic wavej 
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We consider electron waveguides (quantum wires) in graphene created by suitable inhomogeneous 
magnetic fields. The properties of uni-directional snake states are discussed. For a certain magnetic 
field profile, two spatially separated counter-propagating snake states are formed, leading to con- 
ductance quantization insensitive to backscattering by impurities or irregularities of the magnetic 
field. 
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The physics of monolayer graphene devices has re- 
cently attracted a great deal of attention [l], 0] . ^From 
a fundamental perspective, one can hope to relate ex- 
perimental observations to the mathematical properties 
of two-dimensional massless Dirac-Weyl quasiparticles. 
The pseudo-relativistic dispersion relation with Fermi ve- 
locity vp w 10 6 m/sec is intimately connected to the 
sublattice structure: the basis of the graphene honey- 
comb lattice contains two carbon atoms, giving rise to 
an isospin degree of freedom. Graphene has also been 
suggested as new material system for device applica- 
tions 0). In this paper, we pose (and affirmatively an- 
swer) the question whether quantum wires with quan- 
tized conductance can be formed in graphene. Such 
electron waveguides are indispensable parts of any con- 
ceivable all-graphcnc device. In lithographically formed 
graphene 'ribbons', the electronic bandstructure is theo- 
retically expected to very sensitively depend on the width 
and on details of the boundary On top of that, disor- 
der and structural inhomogeneity are substantial in real 
graphene For narrow graphene ribbons or electro- 
statically formed graphene wires Q, conventional con- 
ductance quantization thus seems unlikely [||. This ex- 
pectation is in accordance with recent experiments 0. 

Contrary to such pessimism, we here demonstrate that 
by designing a suitable inhomogeneous magnetic field, a 
magnetic waveguide can be built that indeed allows for 
the perfectly quantized two-terminal conductance 4e 2 //i 
(including spin and valley degeneracy) even when dis- 
order is present. The disorder insensitivity is based 
on a spatial separation of the left- and right-moving 
'snake' states found under the model geometry shown 
in Fig. [Ha). This is reminiscent of the edge states en- 
countered in the integer quantum Hall regime Q, but 
here refers to a completely different microscopic picture. 
Such double-snake states develop in the regime B > but 
B' < 0, while an individual snake state is uni-directional 
and already found in the setup of Fig [Tfb ). Magnetic 
barrier technology is well developed [9L floL [hp ] and its 
application to graphene samples appears to pose no fun- 
damental problems [l^ . In fact, snake states were exper- 
imentally studied in other materials [l3| , mainly mo- 
tivated by the quest for electrical rectification. On the 
theory side, for Schrodinger fcrmions, the magnetic field 
profile in Fig. [Ha) (but only for B' = 0) was discussed 



(homogeneous along y direction), (a) Case <j = 1. For B' < 0, 
counter-propagating pairs of snake states are possible, (b) 
Case a = — 1, with uni-directional propagating snake states. 



in Ref. [14[, and asymmetric cases as in Fig. [IJb) were 
studied by a number of authors [l5| . For the Dirac-Weyl 
quasiparticles encountered in graphene, however, such 
calculations were not reported. Inhomogeneous mag netic 
fields in graphene were discussed by several of us [16[ , and 
we employ that framework in our proposal of magnetic 
waveguides in graphene. 

For a static orbital magnetic field with perpendicu- 
lar component B(x,y), the time-independent Dirac-Weyl 
equation for the quasiparticle isospinor ^{x, y) at energy 
E = vpe reads (we put h = 1) 



-*V- 



(1) 



where following Ref. [16|, we focus on a single K point 
(valley). The Pauli matrices a a with cr = (01,172) act 
in sublattice space, and B(x,y)e z = rotA(a;,y). The 
field profiles considered in Fig. [1] are independent of the 
longitudinal transport direction y and constant within 
each of the three regions, 



B(x) 




(2) 



where a = ±1 gives the relative sign of the magnetic field 
on the two sides |a;| > d. We mention in passing that we 
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have also studied the power-law form B(x) oc x m (with 
m = 1, 2, 3) to make sure that the steps in Eq. do not 
cause unphysical artefacts. Indeed the same qualitative 
features as reported below for the profile |(3J) were found 
from such calculations, which can also benefit from the 
semi-classical approximation. A convenient gauge for the 
vector potential, A = A[x)e y with B(x) = d x A(x), is (for 
<r = +l) 



Bx+{B- B')d, x < -d, 
A{x) = { B'x, \x\ < d, 

Bx — (B — B')d, x > d. 



(3) 



Due to translation invariance in the y-direction, we can 
parametrize solutions fy(x,y) = ip(x)e v by the con- 
served longitudinal momentum k. /From Eq. (JTJ, for 
the spinor component u in ip(x) = (u,v) T , we obtain 



-B 



(x) - (k + ~ c A { x )) + 



0. 



(4) 



For e 7^ 0, v — j^[d x — k — ^A(x)]u then gives the other 
component. To obtain the bandstructure, we first de- 
termine the general solution in each of the three regions 
separately. Matching conditions follow from the continu- 
ity of the wavefunction at x = ^fd and will be shown to 
give an energy quantization condition. 

For x < —d, the constant magnetic field B implies the 
lengthscale l B = y/c/e\B\. We may then explicitly solve 
Eq. (|4|) in terms of parabolic cylinder functions D p (q) 
[la ]. With the auxiliary variables 

q = V2[(x + d)/l B +sgn{B)kl B }, p={el B ) 2 /2-l, (5) 
and complex coefficients a± , the solution reads 



/ D p (±q) 



I D p+1 (±q) 



(6) 
(7) 



Similarly, the cigenfunction for x > d can be expressed 
with coefficients c±, and replacing d — * — d in Eq. l[5]). 
Finally, for B' ^ 0, the region |x| < d again admits such 
a representation with coefficients b± and d — ► in Eq. ([5]). 
For B' = 0, a plane-wave solution applies instead, 



i 



+ik 



a ±ik ± (x+d) 



(8) 



where k± 
is taken as k± 



fc 2 . 



For 



< |fc|, the square root 
Without loss of gen- 



crality we now put B > 0. Normalizability then im- 
plies a+ = c_ a = and we are left with four complex 
coefficients, one of which is fixed by the normalization 
condition. The two matching conditions (at x = ^fd) 
for the 2-spinor ip( x ) then give 4 equations for 3 un- 
knowns, which generates the sought condition for the 
energy bands e n (k). 
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FIG. 2: (a) Spectrum of the magnetic waveguide with 
a = +1, d = Ib and B' — 0. Energies (momenta) are given 
in units of Vf/Ib (Ib 1 )- Only the few lowest electron- like 
(e n (k) > 0) states are shown, (b) Current profile jik F (x) in 
units of vf/Ib, see Eq. with x in units of Ib- The plot 

is for n = 1 and ds = 1, leading to kplB — ±0.7. The two 
counter-propagating states are centered near the middle of 
the waveguide. 



For the symmetric setup a = +1 with B' = 0, some 
algebra yields the energy quantization condition 

W~ 1 (u 2 Vi — Z 2 UiV2) + w(z 2 U 2 Vi — U1V2) (9) 

+ (z 2 - l)(uiu 2 - viv 2 ) = 0, 

which for given k generates an equation for e since fcj_ = 
fcj_(e, k). Here we used the notation 



ui,2 = D p {TV2kl B ), 
/2 

«i, 2 = ±-^—D p+1 (Ty/2kl B ) 
i\e\l B 

w = (k± + ik)/\e\, ; 



(10) 



2ik± d 



Equation ((9]) must then be solved numerically, and leads 
to the energy bands e n (k) shown in Fig. [Ha). For large 
|fc|, the eigenvalues approach the well-known relativistic 
Landau levels at el B = sgn(n)^/2|n| 0], including a zero- 
energy solution (not shown in Fig. 
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FIG. 3: Same as Fig. Ha) but for a = -1, cf. Fig. \V[b). 



To illuminate the current-carrying states, we plot in 
Fig. [Hb) the transverse profile of the particle current, 



jnk (x) = v F {ijj* k (x) ) a 2 ipnk(x), 



(11) 



where tp n k(x) is the transverse eigenspinor to energy 
e n (k). Generalizing the standard argument, see Ap- 
pendix E in Ref. [131, to the case of Dirac-Weyl quasi- 
particles, one can show that 



v n {k) = / dxj nk {x) = d k e n (k). 



(12) 



We stress that Eq. (|12|) is a nontrivial result for Dirac 
fermions. It holds for any magnetic field profile with 
B(x,y) = B(x). This fact leads to the usual can- 
cellation of carrier velocity v n (k) and density of states 
(27r|<9fc£ n (fc)|) -1 , and thus the two-terminal conductance 
will be 4e 2 /ft. (assuming perfect contacts to reservoirs). 
However, as seen in Fig. (2fb), right- and left-moving 
states occupy the same spatial region and are therefore 
susceptible to backscattering perturbations, e.g. due to 
impurities, charge inhomogeneities, or fluctuations in the 
magnetic field. In practice, quantized conductance is 
thus not expected for a waveguide with B' = 0. 

Next we consider the asymmetric case with a = — 1 
but still B' = 0, see Fig. HJb). From the analogy to 
Schrodingcr fermions, one expects to find special uni- 
directional snake states [l5j . On a semi-classical level, 
the uni-directionality can be understood by noting that 
cyclotron orbits have a different winding sense for x < — d 
and x > d. The propagating snake state follows by 
combining half an orbit from each side and a linear 
trajectory in the central region. The energy quantiza- 
tion condition takes again the form © after replacing 

D p+1 (-V2kl B ) and v 2 = j^{p + l)D p (-y/2kl B ) 



"2 



in Eq. (fTU|) . Numerical solution yields the spectrum de- 
picted in Fig. [3l First, we notice a strong asymmetry 
in the energy bands e n {k), just as in the Schrodingcr 
case [ll|. F° r k < a linear dispersion relation is ob- 
served, corresponding to snake states propagating with 
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FIG. 4: (a) Same as Fig.Ha) but for B' = —B, cf. Fig. [Ha). 
The lower pair of e„(fc) curves has an avoided level cross- 
ing (not visible on this scale). The current profile (b) at 
ds = 0.83 (corresponding to fc^'s — ±1) shows that the 
two counter-propagating snake states are spatially separated 
already for d = Ib- 



n- independent velocity \v n \ = vf at sufficiently negative 
k, see Eq. ([T2"j). The equality of snake velocity and Fermi 
velocity for |fce£| 3> 1 also follows from a simple semi- 
classical estimate. Second, the levels merge pairwise at 
large positive k to form the relativistic Landau levels, 
except for the lowest band in Fig. [3] which merges with 
the highest negativc-e band (not shown) to approach the 
zero-energy Landau level. This is a new feature encoun- 
tered only for Dirac fermions and makes this state easily 
identifiable for weakly doped graphene. However, it is 
important to stress that for any finite fc, there can be no 
true zero-energy state for magnetic field configurations 
with cr = — 1. This can be proven on general grounds as 
a consequence of the Aharonov-Casher theorem, which 
in turn follows as a special limit of the celebrated index 
theorem (l8j . 

Interestingly, there is another peculiar subtlety for this 
magnetic field profile. This is seen by computing the 
equilibrium average of the current using Eqs. (jlip and 
(|12|) . which predicts a nonzero result. In fact, the equilib- 
rium current formally diverges and is only limited by the 
bandwidth of the model. To interpret this non-sensical 
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result we note that in the absence of boundaries, the 
snake state propagates in just one direction and thus 
produces an unbalanced current Bow. The conundrum is 
resolved when including boundary contributions to the 
current, which are inevitably present in any real sam- 
ple. In fact, the dispersion relation in Fig. [3] ultimately 
bends upwards for k — ► oo in the presence of a bound- 
ary located at Xf, d. The counter-propagating edge 
state at this boundary will then balance the total cur- 
rent [HI]. We have explicitly checked that this scenario 
holds true for the case of a zig-zag edge, where a sim- 
ple boundary condition on the spinor at x — Xb can be 
used 0]. In analogy to quantum Hall edge states Q, 
however, it should be possible to experimentally probe 
the locally unbalanced current carried by the snake state 
using time-resolved transport measurements fl9j or scan- 
ning tunneling spectroscopy. 

We now go back to the symmetric setup a = +1 but 
take B' < 0. Such a field configuration can be gener- 
ated by depositing two ferromagnetic layers on top of a 
graphene sheet covered by a thin insulating layer [l(| E2] ■ 
In that case one finds two counter-propagating snake 
states, and no boundary contributions are required to get 
zero total current in equilibrium. While for B' = 0, no 
snake states exist, they do appear once B' < 0. By gen- 
eralizing Eq. ([9]), numerical solution of the correspond- 
ing energy quantization condition leads to the results in 
Fig. BJa). Qualitatively, the spectrum consists of snake 
states (with approximately linear dispersion) and Landau 
level states (dispersionless) , with avoided level crossings 
between successive eigenenergies e„(fc). If the Fermi level 



intersects only the lowest band shown in Fig. |4ja), the 
quantized conductance 4e 2 //i follows directly from the 
Kubo formula. The current-carrying states at ±fcp are 
counter-propagating snake states which are spatially sep- 
arated and centered near x = ±d, see Fig. QJb). Due to 
this spatial separation, weak disorder effects or irregu- 
larities in the magnetic field will not be able to induce 
backscattcring processes between these states as long as 
d ^ Ib- In particular, snake states behave identically 
for both K valleys, and thus even inter- valley scattering 
processes are not expected to mix counter-propagating 
states. The conductance quantization in such a setup 
should therefore be observable and very precise. 

To conclude, we have analyzed the properties of elec- 
tron waveguides in graphene, produced by suitable inho- 
mogencous magnetic field profiles. Under the setup in 
Fig. QJa) with B' < 0, we predict robust and highly ac- 
curate conductance quantization in units of 4e 2 /h. This 
system is clearly of interest in the context of interacting 
ID quantum wire physics, as the electron-electron inter- 
action can lead to qualitatively new features. We hope 
that our work motivates experimental and further theo- 
retical studies. 

We thank A. Altland, L. Erdos, T. Hcinzcl and J. Smet 
for discussions. T. K. G. is supported by the A. v. Hum- 
boldt foundation. R. E. is supported by the DFG (SFB 
Transregio 12), and by the ESF network INSTANS. 

Note added: During the preparation of this 
manuscript, a preprint appeared [201 ] where some of our 
results for a = — 1 were also reported. 
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